18 97MNRAS. .57. .342B 


342 


Prof. Brown , Secular Accelerations 


lvii, 5, 


On the Theoretical Values of the Secular Accelerations in the 
Lunar Theory. By Ernest W. Brown. 

(Abstract.) 

The accurate determination of certain inequalities in the 
Moon’s motion has made it possible to find the theoretical values 
of the secular accelerations with a high degree of accuracy. The 
method of determination consists in using the theorem of Professor 
Newcomb for finding the secular accelerations in combination 
with that of the author connecting the mean motions of the perigee 
and node with the constant part of the lunar parallax. The 
resulting values, with the greatest possible amounts by which 
they are uncertain (assuming Newcomb’s value for 8 e' to be cor¬ 
rect), are : 

Secular accelerations per century of the 
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i. The theoretical determination of the secular accelerations, 
formerly one of the most difficult and tedious problems in the 
lunar theory, has been rendered quite simple by the researches of 
Professor Newcomb.* Until the publication of these, two years 
ago, the only investigation with any approach to high accuracy 
was that of Delaunay, who carried the approximations so far as 
to include a great number of terms. The enormous labour of 
finding these, however, was hardly compensated by the result 
obtained, owing to the uncertainty produced by the slow 
convergence of the series when arranged according to powers 
of m. Hansen’s results are entirely numerical, and he used those 
furnished by observation in the construction of his tables. In 
Newcomb’s memoir, just referred to, the values are deduced from 
Delaunay’s algebraical expressions as obtained in the Theorie de 
la Lime, but these again suffer from slow convergence, and errors 
must occur, first, from incorrect estimations of the remainders, 
and, secondly, from errors in Delaunay’s theory. Owing to the 
former cause, Professor Newcomb concludes that his value (5 ;/ *6o) 
for the secular acceleration of the mean motion is uncertain by more 
than Jq of the whole. He adds (mem. cit. p. 199) : “ Although all 
the observations hitherto made on the Moon would not enable us to 
detect errors of this magnitude, a more accurate determination 
would be desirable from a theoretical point of view.” It is the 

* Papers published for the use of the Amer. Eph. vol. v. pt. iii. 
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object of this paper to show how such a determination may be 
made, and to find the values correctly to 3^ part of the whole 
in each case. 

This has been rendered possible in three ways : First, by 
Newcomb’s theorem, showing how Delaunay’s final constants 
L, G, H may be calculated from the expressions for the rect¬ 
angular coordinates of the Moon ; secondly, by my theorem 
connecting the constant term in the expression for the parallax 
with L, G, H and the motions of the perigee and node ; * and, 
thirdly, by the new determinations of these mean motions, the results 
of which have been given in my paper “On the Mean Motions of 
the Lunar Perigee and Node.” t These theorems have also been 
employed to test Delaunay’s algebraical expressions for the secular 
accelerations, and the results are given in 4 4 below. In order to 

0 2 7T 

find the numerical values, such expressions as r -—have to be 

tino.e ' 2 

calculated ; it is shown in § 5 how they may be found with very 


fair accuracy when the complete numerical value of 


0 


7r 


de ,%d 


and 


only a few terms of its slowly convergent algebraical expression are 
known. Finally, the numerical results are given, with the 
amounts by which they are uncertain. 


General Formulce for the Secular Accelerations. 


2. Professor Newcomb’s theorem (mem. cit. p. 191) states that 
the variations of n (or a), e , y due to a variation he' of e' may be 
obtained accurately to the first order, with respect to S, by solving 
the equations 


or 


8L = o. 8Gr = o, 2H = o. 

0L ? 0L * - 3l . „ ?L * ,* p 

5 ft + r _ 8 . e 2 + -—5.7-=— 5 .e r2 , &c.; 


B .e 2 


0 . y- 


0 . ( 


and that the secular accelerations of the mean motion, the perigee 
and node, are then given by 


§e= Jsft dt, 87r = j'57r 1 dt, 50 = 1 * 20 , dt 3 

where tt,, 6 1 are the mean motions of the perigee and node. 
Putting L=Cj, G — L—c 2 , H —G=c 3 , the first three equations 
become 

8c x = 0 , 5 c 2 = o, 5 c 3 == o. 

* Proc. London Math. Soc. 1897. 
t Monthly Notices Ii.A.S. 1897. 
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I have shown in the latter of the two papers just referred to 
that the variations of n , e , y, e' also satisfy the equations 


ft, S ,1^+ ft, 8 ■ =1 (E + M) S . 

2 On 


d, 


d, 


©.• 


Cn 5 


3*. +c 8 39, 

* ' 9 -r 


= 2 (E + M) 5 . JL 
2 d • V 2 \ rJ o 


These results immediately follow from those given in my 
former paper in the Monthly Notices by submitting them to a 
variation 8 and putting 8c 1} Sc 2 , 8c 3 zero. As before, E + M =n 2 a 3 , 

and (—j 0 denotes the constant term in the expression of 1 ; the 
\rj r 

constants are defined as in Delaunay’s final expressions for the 
coordinates. The square of the ratio of the solar to the lunar 
parallax is neglected, but its effect is quite insensible. 

The first of these three equations may be used effectively to 
find the secular accelerations. We have 


07 T, _ 0 2 7 Tj 

’cT 


n Jn* 


$n + 


d 2 i r, s , 9 , 0-7T, s . o , 0‘V, * 

5 e + 5rV"d 5 5 • 7-, 

On 0 • e - on O • e- On 0 ■ y- 


with similar results for the other expressions. It can be easily 
shown that c 2 , 8*e 2 contain the factor c 2 in all their terms, and 
c 3 , 8*y 2 the factor y 2 . Also it is known by a theorem of Adams * 

that there are no terms in ( - ] of the form - e 2 f (m.e' 2 ), 

\rj o a 

1 y 2 / Hence, in order to find Zn to the orders e 2 , y 2 , and 


a 


all powers of e' 2 , we can neglect 8 . c 2 , 8 . y 2 in the first equation, 
so that it becomes 

ft, (t—’ 5 *+a% L r i 5 - e ' 2 ) + c » (1^ 5?i + fTX L T s - 

\On 2 OnO.e ' 2 / \on 2 OnO.e 2 


= ? (E + M) in + W* ] , 

2 L dn 2 9 n d.e“ J 

putfor (i) t . 


3 2 p« 


where P 0 has been 

If we neglect e 2 , y 2 , it is seen that the principal part of Zn is 
given by 

H2D r\ 2D 

« 5/ 2 ....... (I) 


3 2 P„ fa ^ 0 2 P. 

0?i 2 3 « 0 .e' 2 


* Monthly Notices f vol. xxxviii. pp. 460-472. Coll. Works, vol. i. pp. 189-204. 
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When the principal part of 8n has been found from this 
equation, its value may be substituted in the left-hand member 
of the previous equation, and the parts depending on e 2 , y 2 can 
then be evaluated. 

For the motion of the perigee we have 

fcr.aA'>-S. e'°- +^Sn +5 . $ . y\ 

o • e ' z d n c • f a • 7 2 


If we neglect e 2 , y 2 , it will be sufficient to use the first two terms 
only. The parts of chief importance in all cases arise from 
07 Tj /9. e ' 2 , and their accuracy depends mainly on the accuracy with 
which the terms having the factor e’~ in - l are known. If, how¬ 
ever, the parts in ffirj, depending on e~. y-, be required, it will be 
simplest to find 8. e 2 , 8 . y 2 from the equations 8c 2 =o, 8c 3 =o, 
respectively. The principal parts of 8. e~, 8 . y 2 are by this method 
found in the forms 


^8 . e 2 — e-f '(nCfin, n$ . y- — y-f 


and can therefore be found, since the value of 8 n is already 
known. 

In the present state of the lunar theory nothing is to be 
gained by calculating these higher portions in the case of the 
perigee, for the complete numerical values of the terms in 7r ]? 
wffiich have the factors eV 2 , yV 2 have not yet been obtained, 
and Adams’ value of P 0 furnishes only the terms up to the orders 
e‘V 2 m 3 , yV 2 ra 3 . Hence we must verify Delaunay’s values to 
this order and estimate for the remainders. Fortunately, owing 
to the smallness of the quantities involved, the possible errors of 
estimation are very small and can be easily seen. 

Exactly similar results hold for 8 ( 9 ,. 

3. The equation (1) explains partly the slowness of the con¬ 
vergence of the expression which represents the principal part of 
8w when expanded in powers of m. All inequalities in the 
lunar theory which contain the factor e' appear to be more or less 
slowly convergent. When such expressions as e' 2 f(m) are 
differentiated with respect to n , the earlier terms will converge 
more slowly still, but the convergence of the later terms will not 
be much affected. The calculations further show that slow 
convergence in the higher terms will arise partly from the 

3 2 P 

expansion of the inverse of in powers of m ; a remedy for this 

is evident. The convergence of all the results is, as usual, 
improved by using m { =-n f /(n—n') instead of m—n’jn as a 
parameter. 

4. A Test of Delaunay’s Expressions for the Secular 
Accelerations .—These are given in C.E. vol. xliv. p. 209, in terms of 
Plana’s constants, and, with an increased number of terms, in C.R. 
vol. lxxiv. p. 174, in terms of his own constants. The constant part 
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of i j.r has been obtained by A dams (Coll. Works , vol. i. p. 203), and 
his value was used in the formulae of § 2. The values of c 2 , c 3 are 
given by Newcomb (mem. cit. p. 184), and the values of tt u 0 x 
are those of Delaunay ( C.B. vol. lxxiv. p. 19) with the corrections 
noted in my paper in the Monthly Notices already referred to 
(§ 1). The last powers of on tested in each ease are as follows :— 
The value of Se was tested as far as m 7 , y 2 m 5 , ehW, e'hrC, and 
the values of Sw, SO as far as m 5 , y 2 m 3 , e 2 m 3 , e' 2 m 2 . The only 
differences found were the coefficients of the terms involving y 2 m 3 


in Sw and e 2 m 3 in SO, for which 


' j ——7^ were found instead 
4 16 


of —• . The terms involving a 2 j a ' 2 and higher powers 

4 16 

and products of e 2 , y 1 , e' 2 do not produce additions so great as 
o"*oi, and they will therefore not be considered. 

5. Method of combining algebraical and numerical results to 
obtain the derivatives with regard to n.—The problem may be 
stated as follows : Given the first i terms and also the complete 
numerical value of a series whose terms are diminishing in abso¬ 
lute value, to find a close approximation to the value of its 
derivative. The following is a practical method of solution :— 
Let 


f[m) = a 0 + a,m + a^mr + . . . , 


where a 0 , a 1? . . are numerical coefficients, and put 
fi( m ) = a o + a \ m + a 2 m ~ + * • • 


Taking the derivatives with regard to on, we can write the 
result 


_i_ j 

agn 1 + --- a i+1 m i+l 


+ 



If i be not too small (e.g. not less than 5), and the terms of the 
series are not diminishing too slowly (e.g. with a ratio not much 
greater than ^), the error committed by putting 


mf(m) = - mffn)] 

will be small compared with the true value of mf{gn). We shall 
generally obtain greater accuracy by using i-f- a proper fraction 
instead of i as the factor, when the series appears to be diminish¬ 
ing with fair regularity. 

For example, the part of 774 which depends on e ' 2 when the 
parameter m Y is used, is 


nc n ( -m x 2 + m* + ^ 22 43 m \ + 94 2 777 w s\ 

\8 32 1 256 1024 ) 


Here i= 6 , and reducing the portion in brackets to numbers, we 
obtain 

+ *00735+ *012444 *00705 +-00318= + *03002. 
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I find 
Hence 


/K)= +-032063. 


A m i ) \) W : + *00204. 


Applying the formula, and using 6^ as the multiplying factor 
of this difference, we find 


f(m } ) = + -0147 + -0373 + -0282 + *o 159 + '0133 - + -1094. 


The result might possibly be as much as four or five units 
wrong in the last place of decimals, but in all the cases treated 
below, the possible error does not affect the results to the degree 
of accuracy required. To complete the determination, we have, 
for this portion of 

I s -(* +" 

on 

- + *0321 --1182=- *oS6, 


which is not more than one unit wrong in the third place of 
decimals. 


The Numerical Values. 

6 . The following are the numerical data required :— 

In the values of c 2 , c 2f we only require the terms which are 
of the forms na 2 e 2 f(m) i na 2 y 2 f(m) i respectively. These may be 
accurately obtained by Newcomb ? s method from the expressions 
in terms of rectangular coordinates, c 3 from the results in a 
paper by Mr. P. H. Cowell,* and c 2 from those in a paper by the 
writer.! We find 


<7 2 = _ -47 4 o35^ 2 A c 3 — — 2002 S 7 nct - y 2 . 

The values of 7r l5 are as follows : i 

ttj = n [ + *0085726 + -032063c' 2 ], 

^^[ — •0039992 —-o°5224e /2 J. 


Proceeding as in § 5, we then find 

—.= —-01502 — *o86e' 2 , 
c)n 


de, 

on 


= + -003708 + -00319c' 2 , 


d-j r. o 

n = — *° 4 ^ 9 j 

On* 


3 % . 

n = — 


00661. 


* Amer. Jour . Math. vol. xviii. p. 113. 
t Ibid . vol. xv. p. 261. 

j See the paper “ On the Mean Motions of the Lunar Perigee and Node.*’ 

D D 
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The numbers contain powers of m only ; the omitted portions, 
depending on e 2 , etc., are not required. 

Further, from Adams’ algebraical results and my numerical 
calculations : 


aP Q = 1*00091 + *ooii75oe' 2 + *000640' 4 , 

3p 

an = + •66551 — *ooo9865^ 3 + *002 6 e'\ 

on 


an 2 



— •21975 — *000131 $e’K 


7. Pursuing the operations of § 2 above, we obtain 

hi =[-• 004489 + ’024e ,2 + • 123 5e 2 - *oi967 2 }m5 . e’\ 
d 7 r t = + *03213^5 , e' 2 , 

—*005241^5 . e f % 


Put now, with Delaunay, 

(n 5 e r 2 dt= —1270"!% 


where T is the time reckoned in centuries of Julian years, and 
include the parts of Stt 15 80 1 which depend on e 2 , y 2 , e n from his 
expressions (see § 4 above), estimating for the remainders in each 
case. The maximum possible error due to the latter is appended, 
with the resulting values in the table : 


Secular Accelerations. 


Parts of form. In Bn. 

f(m) +570 

—0*01 

e 2 f{m) +0*47 

Tf(m) -0*05 

Complete values 4 - 6* 11 ± o"*02 


In Sttj. 

In B9 t . 

n 

it 

— 40*8 

+ 6-66 

+ 0*0 

— 000 

+ 0*2 

±0*14 

+ 0*4 

— 002 


-40*2±0' / *i +6*78 ±o"*o2 


These results agree almost exactly with those of Delaunay. 
With the latest values for the masses of the planets, Professor 
Newcomb finds (mem. cit. p. 199), 

|a' 2 = -1229"!% 


Hence the definitive theoretical values of the secular accelera¬ 
tions per century are, for 


The Mean Motion 
The Perigee 
The Node 


n n 

+ 5 91 ± 0*02; 
— 38*9 ±o*l; 
+ 6*56 ±002. 
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Laying aside the question of the accuracy of the value of 
j%8e' 2 , the above theoretical values for the secular acceleration 
are correct to 3!^ part of the whole in each case. 

Haverford College , Pa.: 

1897 January 11. 


On a Photographic Transit Circle. By H. H. Turner, M.A,, 
B.Sc., Savilian Professor. 

1. The desirability of introducing photography into meridian 
observations, as a method of eliminating personal and other errors, 
has been recognised by many astronomers, and will probably 
be admitted by all, if this introduction can be accomplished with¬ 
out endangering the simplicity and stability of the present form 
of instrument. The following solution of this problem, which 
presents some advantages not realised in other forms, is offered 
for consideration, 

2. It may be well in the first instance to enumerate the chief 
features of the present instrument—its advantages and disadvan¬ 
tages. 

Features : 

(a) Movable visual telescope with two micrometers 
at right angles. 

(b) Two carefully made pivots, allowing one degree of 
freedom only. 

(c) Graduated circle and microscopes. 

(d) Chronograph. 

All these, or their equivalents, may thus be allowed in any new 
arrangement. 

Advantages : 

(e) Stability, owing to there being only one degree of 
freedom, excepting in so far as the micrometers are used. 

(/) Simplicity. 

The advantages may generally be brought under one or other 
of these heads. It is unnecessary to specify them in detail unless 
we are going to propose any rearrangement of the instrument 
which will sacrifice any of them, and that is not my present 
purpose. It is more important to consider carefully how we may 
remove the 

Disadvantages of the Present Instrument, viz. 1 

(g) Personal equation of all kinds. 

(h) Flexure and irregular refraction. 

( i) Limitation to comparatively bright stars. 

D D 2 
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